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Abstract. Critical behaviour of two systems, subjected to the turbulent mixing, 
is studied by means of the field theoretic renormalization group. The first system, 
described by the equilibrium model A, corresponds to relaxational dynamics of a 
non-conserved order parameter. The second one is the strongly non-equilibrium 
reaction-diffusion system, known as Gribov process and equivalent to the Reggeon field 
theory. The turbulent mixing is modelled by the Kazantsev-Kraichnan "rapid-change" 
ensemble: time-decorrelated Gaussian velocity field with the power-like spectrum 
oc k~ d ~t. Effects of compressibility of the fluid are studied. It is shown that, depending 
on the relation between the exponent £ and the spatial dimension d, the both systems 
exhibit four different types of critical behaviour, associated with four possible fixed 
points of the renormalization group equations. Three fixed points correspond to known 
regimes: Gaussian fixed point, original model without mixing and passively advected 
scalar field. The most interesting fourth point corresponds to a new type of critical 
behaviour, in which the nonlinearity and turbulent mixing are both relevant, and 
the critical exponents depend on d, £ and the degree of compressibility. The critical 
exponents and regions of stability for all the regimes are calculated in the leading 
order of the double expansion in two parameters £ and e = 4 — d. For the both models, 
compressibility enhances the role of the nonlinear terms in the dynamical equations: 
the region in the e— £ plane, where the new nontrivial regime is stable, is getting much 
wider as the degree of compressibility increases. For the incompressible fluid, the 
most realistic values d = 3 and £ = 4/3 (Kolmogorov turbulence) lie in the region of 
stability of the passive scalar regime. If the compressibility becomes strong enough, 
the crossover in the critical behaviour occurs, and these values of d and £ fall into 
the region of stability of the new regime, where the advection and the nonlinearity are 
both important. In its turn, turbulent transfer becomes more efficient due to combined 
effects of the mixing and the nonlinear terms. 
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1. Introduction 

Numerous systems of very different physical nature exhibit interesting singular 
behaviour in the vicinity of their critical points (second-order or continuous phase 
transitions). Their correlation functions reveal self-similar behaviour with rather 
universal critical dimensions: they depend only on few global characteristics of the 
system (like symmetry or space dimensionality). Quantitative description of critical 
behaviour is provided by the renormalization group (RG). In the RG approach, possible 
types of critical regimes (universality classes) are associated with infrared (IR) attractive 
fixed points of renormalizable field theoretic models. Most typical equilibrium phase 
transitions belong to the universality class of the On-symmetric ip 4 model of an n- 
component scalar order parameter. Universal characteristics of the critical behaviour 
depend only on n and the space dimensionality d and can be calculated in the form 
of the expansion in e = 4 — d or within other systematic perturbation schemes; see 
the monographs [U [2] and the literature cited therein. Equilibrium dynamical critical 
behaviour is much richer and less universal; it is described by variety of models; see the 
reviews [3] and chap. 5 of book [2]. 

Over the past few decades, much attention has been attracted by the spreading 
processes and corresponding nonequilibrium phase transitions; see e.g. [U |5] and the 
literature cited therein. Spreading processes are encountered in physical, chemical, 
biological, ecological and sociological systems: autocatalytic reactions, percolation in 
porous media, forest fires, epidemic diseases and so on. Depending on the conditions, 
the spreading of an "agent" (e.g. infectious disease) can either continue over the whole 
population, or terminate after some time. In the first case, the system evolves to a 
stationary (but not thermally equilibrium) active state, in which the density of the 
agent remains a fluctuating quantity. In the second case, the system is "trapped" in 
an absorbing (inactive) state where all the fluctuations cease completely. Transition 
between these phases are continuous; they are especially interesting as examples of 
nonequilibrium critical behaviour |4j. The most popular model that describes the 
spreading of an agent (e.g., infectious disease) in a reaction-diffusion system belongs 
to the universality class of the directed bond percolation process, also known as simple 
epidemic process or Gribov process pE]. In the field theoretic formulation that model 
is equivalent to the Reggeon field theory and was also studied within the RG approach 
and the e expansion [5]. 

However, it has long been realized that the behaviour of a real system near 
its critical point is extremely sensitive to external disturbances, geometry of the 
experimental setup, gravity, presence of impurities and so on; see e.g. the monograph 
[6] for the general discussion and the references. "Ideal" critical behaviour of an infinite 
stationary system can be obscured by finite-size effects, finite time of evolution (ageing) 
and so on. What is more, some disturbances (randomly distributed impurities or 
turbulent mixing) can produce completely new types of critical behaviour with rich 
and rather exotic properties, like e.g. expansion in \Je rather than in e E]. 
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Those issues become especially important for the nonequilibrium phase transitions, 
because the ideal conditions of a "pure" stationary critical state can hardly be achieved 
in real chemical or biological systems, and the effects of various disturbances can never 
be completely excluded. In particular, intrinsic turbulence effects cannot be avoided 
in chemical catalytic reactions or forest fires. One can also speculate that atmospheric 
turbulence can play important role for the spreading of an infectious disease by flying 
insects or birds. 

Investigation of the effects of various kinds of deterministic or chaotic flows (laminar 
shear flows, turbulent convection and so on) on the behaviour of the critical fluids (like 
binary liquid mixtures near their consolution points) has shown that the flow can destroy 
the usual critical behaviour, typical of the if) A or Reggeon models. It can change to the 
mean-field behaviour or, under some conditions, to a more complex behaviour described 
by new nonequilibrium universality classes [8] |15j. 

In this paper we study effects of turbulent mixing on the critical behaviour of two 
systems near their critical points, paying special attention to compressibility of the fluid. 
To describe the phase transitions, we use two paradigmatic models of dynamical critical 
behaviour. The first one is the equilibrium model A, which describes purely relaxational 
dynamics of a non-conserved order parameter (e.g. in a nematic liquid crystal); see the 
reviews [3] and chap. 5 of book |2]. The second one, referred to as Gribov model in 
what follows, is the aforementioned model of a nonequilibrium reaction-diffusion system 
HE]. 

To model the turbulent mixing, we employ the famous Kazantsev-Kraichnan 
"rapid-change" ensemble: time-decorrelated Gaussian velocity field with the pair 
velocity function of the form (vv) oc 8(t — if) k~ d ~^, where k is the wave number and 
< £ < 2 is a free parameter with the most realistic ( "Kolmogorov" ) value £ = 4/3. 
This model has attracted enormous attention recently because of the insight it offers 
into the origin of intermittency and anomalous scaling in the fully developed turbulence; 
see the review paper [16] and references therein. The RG approach to that problem 
is reviewed in [17J. In the context of our study it is especially important that the 
Kazantsev-Kraichnan ensemble allows one to easily model compressibility of the fluid, 
which appears rather difficult if the velocity is modelled by dynamical equations; see 
e.g. discussion in [T8] . 

The outline of the paper and the main results are the following. In section [2] we 
present the detailed description of the models and their field theoretic formulation. In 
section [3] we analyze the ultraviolet (UV) divergences, relaying upon the power counting 
and additional symmetry considerations. We show that the both models, after proper 
extension, are multiplicatively renormalizable. Thus we can derive the RG equations and 
introduce the RG functions (/? functions and anomalous dimensions 7) in the standard 
fashion; see section HI 

In section [5] we show that, depending on the relation between the spatial dimension 
d and the exponent £ in the velocity correlator, the both models exhibit four different 
types of critical behaviour, associated with four fixed points of the corresponding 
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RG equations. Three fixed points correspond to known regimes: ordinary diffusion 
(Gaussian or free field theory), non-interacting scalar field passively advected by the flow 
(the nonlinearity in the order parameter in the original dynamical equations appears 
unimportant), and the original critical behaviour of a model without mixing (that is, 
model A or Gribov model). The most interesting fourth point corresponds to a new full- 
scale nonequilibrium universality class, in which both the nonlinearity and turbulent 
mixing are relevant. The corresponding critical exponents can be systematically 
calculated as double expansions in two parameters: £ and e = 4 — d. 

The practical calculation of the renormalization constants, RG functions, regions 
of stability and critical dimensions was accomplished to the leading order (one-loop 
approximation); some of the results, however, are exact (valid to all orders of the 
double e— £ expansion). In general, the critical exponents and the boundaries between 
the regions of stability of critical regimes in the e— £ plane depend on the degree of 
compressibility. 

The results obtained are discussed in section [7J The main qualitative conclusion 
of our analysis is that, for the both models, compressibility enhances the role of the 
nonlinear terms in the dynamical equations. The region in the £-£ plane, where the new 
nontrivial regime is stable (the corresponding fixed point is attractive), is getting much 
wider as the degree of compressibility increases. For the case of incompressible fluid, 
the most realistic values al = 3 and £ = 4/3 lie in the region of stability of the passive 
scalar regime, that is, the nonlinear terms are irrelevant for the critical behaviour. If the 
compressibility becomes strong enough, the crossover in the critical behaviour occurs, 
and these values of d and £ fall into the region of stability of the new regime, where the 
advection and the nonlinearity are both important. 

These general considerations are illustrated by the example of a cloud of particles, 
randomly walking in a nearly critical turbulent medium. For a weakly compressible 
fluid, spreading of such a cloud is determined solely by the turbulent transfer. If the 
compressibility becomes strong enough, the new scaling regime comes into play and 
the transfer becomes faster, due to combined effects of the mixing and the nonlinear 
terms. The corresponding exponent depends explicitly on the degree of compressibility 
and differs from the well-known "1/2 law" for ordinary diffusion and from Richardson's 
"4/3 law" for turbulent transfer. 

2. Description of the models. Field theoretic formulation 

In the Langevin formulation the models are defined by stochastic differential equations 
for the order parameter ip = ip(t,x): 



where d t = d/dt, d 2 is the Laplace operator, Ao > is the kinematic (diffusion) 
coefficient and tq oc (T — T c ) is the deviation of the temperature (or its analog) from 
the critical value. The nonlinearity has the form V(ip) = Uoip 3 /3\ for the model A and 



d t ip = A {(-r + d 2 )4> - V(ip)} + C = 0, 
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V(ip) = g ip 2 /2 for the Gribov process; g Q and u > being the coupling constants. The 
Gaussian random noise £ = x) with the zero mean is specified by the pair correlation 
function: 

(C(t,x)C(t',x')) =2A 5(t-i')5 (d) (x-x') (2.2) 

for the model A and 

(C(t, x)C(t', x')> = SoAo V(* , x) - t')5^ (x - x') (2.3) 

for the Gribov process; d is the dimension of the x space. The factor ip in front of the 
correlator ( 12. 31) guarantees that in the absorbing state the fluctuations cease entirely, 
while the factor 2Ao in (12.21) ensures the correspondence to the static if) 4 model. Here and 
below, the bare (unrenormalized) parameters are marked by the subscript "o." Their 
renormalized analogs (without the subscript) will appear later on. 

The stochastic problems (12. ip . (12. 2p . (12. 3p can be reformulated as field theoretic 
models of the doubled set of fields $ = {ip, ip^} with action functional 

Sty, tf) = ^ (-d t + X d 2 - A r ) V + A (^) 2 - u ^if) 3 /3\ (2.4) 
for the model A and 

Sty, tf) = 1>H-b + ^od 2 - \ r )i) + ^° {(^) 2 ifj - ^V 2 } (2.5) 

for the Gribov model. Here, ip^ = V* (^> x ) i s the auxiliary "response field" and the 
integrations over the arguments of the fields are implied, for example 

^d t ip = I dt I rfx^(t,x)^(t,x). 



The field theoretic formulation means that the statistical averages of random quantities 
in the original stochastic problems can be represented as functional integrals over the full 
set of fields with weight exp5($), and can therefore be viewed as the Green functions 
of the field theoretic models with actions (12. 4ft . (12. 5ft . In particular, the linear response 
function of the problems (12. ip . (12. 2p . (12. 3p is given by the Green function 



G= (V> t (t,x)V>(t / ,x')) = J V ^J V ^ ^'XW,*') expS{ip,ip r ) (2.6) 

of the corresponding field theoretic modeljij] 

The model (12. 5p corresponds to the standard Feynman diagrammatic technique 
with the only bare propagator (ipip^o and two triple vertices ~ (ip^) 2 ip, ip^ip 2 . In the 
time-momentum and frequency-momentum representation the propagator has the form 

(W*)a(t, k) = 6{t) exp {-A (A; 2 + r )t} , 

*) = ■ nAr^T y (2-7) 
—iu + A (Ar + T ) 



| The continuous (^-correlated stochastic noise like £ in (|2.3j) . usually referred to as multiplicative 
random noise, requires additional careful definition; see e.g. |19j . Consistent description is provided by 
the so-called Ito time discretization, which is also assumed here. 



Effects of turbulent mixing on critical behaviour in the presence of compressibility 6 

Here 9{. . .) is the Heaviside step function, so that the propagator (12. 7p is retarded^ 
Then from the analysis of the diagrams one can check that the Green functions built 
solely from the field if> or solely from if^ necessarily contain closed circuits of retarded 
propagators (12.71) and therefore vanish identically. 

For the functions (ifv . . . if) ') this fact is a general consequence of the causality, which 
is valid for any stochastic model; see e.g. the discussion in [2j. Then vanishing of the 
functions (if . . . iff) in the model (12.51) can be viewed as a consequence of the symmetry 
with respect to the transformation 

V>(t,x) -^t(_ t) _ x ) 5 ^t( fjX )_ > ^,(_ tj _ x ) j g ^-g Q . (2.8) 

Reflection of the constant go is in fact unimportant because, as can easily be seen, in 
the model (12. 5p with trilinear interaction terms the actual expansion parameter in the 
perturbation theory is g\ rather than g itself; we shall also denote it as u = g\ in what 
follows. 

In addition to (12.71) . diagrammatic technique for the model A involves the 
propagator (ipip)o of the form 

(if)if))o{t,k) = 1 exp{-A (A; 2 + r )|t|}, 

2A 



while the only vertex has the form ~ ip^if) 3 . 

The Galilean invariant coupling with the velocity field v = {t>j(t,x)} for the 
compressible fluid (diVi ^ 0) can be introduced by the replacement 

d t ip -> d t ip + a di(viip) + (a - l)(vidi)ip = V t ip + a (diVi)ip (2.10) 

in (12.11) . Here V< = dt + Vidi is the Lagrangian derivative, ao is an arbitrary parameter 
and di = d/dxi. For the linear advection-diffusion equation, the choice ao = 1 
corresponds to the conserved quantity if) (e.g. density of an impurity) while ao = 
is referred to as "tracer" (concentration of the impurity or the temperature of the fluid; 
in this case, the conserved quantity is if)'). As we shall see below, in the presence of 
nonlinearities in ( 12. ip . it is necessary to simultaneously include the both terms (12.101) in 
order to ensure multiplicative renormalizability. 

In the real problem, the field v(t, x) satisfies the Navier-Stokes equation. We will 
employ the rapid-change model [16], where the velocity obeys a Gaussian distribution 
with zero mean and the correlation function 

(^xK^x'^^-OAiW, r = x-x' (2.11) 

with 

A,(r) = D {P«(k) + aQ^k)} exp(ikr). (2.12) 

S The Ito discretization accepted above means that all the diagrams with self-contracted (ip'ip)o lines 
should be discarded. Vanishing of the propagator (ip^ip^o is the general feature of stochastic dynamical 
models. 



Effects of turbulent mixing on critical behaviour in the presence of compressibility 7 



Here i^j(k) = 5^ — kikj/k 2 and Qij(k) = kikj/k 2 are the transverse and the longitudinal 
projectors, k = |k| is the wave number, D > is an amplitude factor and a > is an 
arbitrary parameter. The case a = corresponds to the incompressible fluid (diVi = 0), 
while the limit a — > oo at fixed aD corresponds to the purely potential velocity field. 
The exponent < £ < 2 is a free parameter which can be viewed as a kind of Holder 
exponent, which measures "roughness" of the velocity field; the "Kolmogorov" value is 
£ = 4/3, while the "Batchelor" limit £ — > 2 corresponds to smooth velocity. The cutoff 
in the integral ( 12.12ft from below at k = m, where m = 1/ C is the reciprocal of the 
integral turbulence scale C, provides IR regularization. Its precise form is unimportant; 
the sharp cutoff is the simplest choice for the practical calculations. 

The action functionals for the full set of fields $ = {tp, ip>, v} become 

S A ($) = <S C ($) + A (^) 2 - « Ao^73! (2.13) 
for the model A and 

Sap) = Sc{*) + ^ {(^ - ^V 2 } (2.14) 

for the Gribov model, where the common part of the actions has the form 

S c m = ^ {-V t + A {d 2 - r ) - a (d iVl )}^ + 5(v). (2.15) 

They are obtained from (12.41) . (12.51) . by the replacement (I2.10p and adding the term 
corresponding to the Gaussian averaging over the field v with the correlator (12. lip . 
d2l2D: 



^JaltJdxJ dx! Vi {t,x)D^{r) Vj {t,x!), (2.16) 



where 



5v = - 

v ; 2 



is the kernel of the inverse linear operation for the function -Dij(r) in ( 12. 12ft . 

In addition to (|2.7|) . the Feynman diagrams for the models (I2.13I) - (I2.15I) involve 
the propagator (vv)o specified by the relations ( 12.111) . ( 12.12ft and the new vertex 

■tfviVd) = {(vidi)ij + aoidiVi)} i>. (2.17) 

In the diagrams it corresponds to the vertex factor 

Vi = -iki - ia <7i, (2.18) 

where ki is the momentum argument of the field ip and is the momentum of t>i. 

The reflection symmetry ipi^ ~^ ~i J , °f the model A is preserved in the full 
model ( 12.131) . while for the full model ( 12.141) the relations ( 12. 8p should be augmented by 
the transformation 

ao->>(l-ao), (2.19) 
as can directly be seen from (12.151) with the aid of integration by parts. 
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For the both full models, the role of the coupling constants (expansion parameters 
in the perturbation theory) is played by the three parameters 

u ~ A 4 ~ d , w = D /X ~ A 5 , w a ~ A ? . (2.20) 

The last relations, following from the dimensionality considerations (more precisely, 
see the next section), define the typical UV momentum scale A. From the relations 
( I2.20p it follows that the interactions i/j'i/j 3 , (if)^) 2 ^ and i/j^i/j 2 become logarithmic (the 
corresponding coupling constant uq becomes dimensionless) at d = 4. Thus for the 
single-charge problems ( 12. 41) . ( 12. 51) . the value d = d c = 4 is the upper critical dimension, 
and the deviation e = 4 — d c plays the part of the formal expansion parameter in the 
RG approach: the critical exponents are nontrivial for e > and are calculated as series 
in e. 

The additional interactions ~ ip^vdip of the full models ( 12.13!) . (12. 14}) become 
logarithmic at £ = 0. The parameter £ is not related to the spatial dimension and 
can be varied independently. However, for the RG analysis of the full problems it is 
important that all the interactions become logarithmic at the same time. Otherwise, 
one of them would be weaker than the others from the RG viewpoint and it would be 
irrelevant in the leading-order IR behaviour. As a result, some of the scaling regimes of 
the full model would be overlooked. In order to study all possible scaling regimes and 
the crossovers between them, we need a genuine three-charge theory, in which all the 
interactions are treated on equal footing. Thus we will treat e and £ as small parameters 
of the same order, e cx £. Instead of the plain e expansion in the single-charge models, 
the coordinates of the fixed points, critical dimensions and other quantities will be 
calculated as double expansions in the e-£ plane around the origin, that is, around 
the point in which all the coupling constants in f l2.20f) become dimensionless. Similar 
situation was encountered earlier in various models of turbulence and complex critical 
behaviour, e.g. P2 Q31 El H51 12B] ■ 

3. Canonical dimensions, UV divergences and the renormalization 

It is well known that the analysis of UV divergences is based on the analysis of canonical 
dimensions ("power counting"); see e.g. [U(2]. Dynamical models of the type (12. 4J) , ( 12. 5 j) 
and (I2.13p -l l2.15p . in contrast to static ones, have two independent scales: the time scale 
T and the length scale L. Thus the canonical dimension of any quantity F (a field or 
a parameter) is characterized by two numbers, the frequency dimension d F and the 
momentum dimension dp, defined such that [F] ~ [T] ~ d F [L] ~ dp . These dimensions are 
found from the obvious normalization conditions 

d\ = -d k x = 1, dt = d£ = 0, d* = d\ = 0, c = -d? = 1 

and from the requirement that each term of the action functional be dimensionless (with 
respect to the momentum and frequency dimensions separately). Then, based on dp and 
d F , one can introduce the total canonical dimension dp = dp + 2d F (in the free theory, 
d t oc d 2 ), which plays in the theory of renormalization of dynamical models the same 
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Table 1. Canonical dimensions of the fields and parameters in the models ([2.130 - 
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A , A T , T 
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it, w, a, ao, a 
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al/2 - 1 


d/2 + 1 


al/2 


-1 


-2 2 
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4-d 







d" F 











1 


1 














oIf 


al/2 - 1 


al/2 + 1 


al/2 


1 


2 


1 


4 - al 








role as the conventional (momentum) dimension does in static problems; see Chap. 5 of 
[2]. The canonical dimensions for the models ( 12 . 1 3 [) — ( ]2 . 1 5 [) are given in table HJ including 
renormalized parameters (without subscript "o"), which will be introduced soon. The 
fields of the model A are marked by the subscript A, while the fields of the Gribov 
model are marked by G. The dimensions of the parameters Ao, To etc are the same for 
the both models. 

As already discussed in the end of the previous section, the full models are 
logarithmic (all the coupling constants are simultaneously dimensionless) at d — 4 and 
£ = 0. Thus the UV divergences in the Green functions manifest themselves as poles in 
e = 4 — d, £ and, in general, their linear combinations. 

The total canonical dimension of an arbitrary 1-irreducible Green function T = 
($ • ■ ■ $)i_i r is given by the relation [2] 

d T = 4 + 2d£ = d + 2 - N^d®, (3.1) 

where N$ = {N^,N^,N V } are the numbers of corresponding fields entering into the 
function T, and the summation over all types of the fields is implied. The total dimension 
dr in logarithmic theory (that is, at e = £ = 0) is the formal index of the UV divergence 
Or = dr| e =5=o- Superficial UV divergences, whose removal requires counterterms, can 
be present only in those functions T for which o~r is a non-negative integer. 
From table Q] and (13.1 ft we find 

5 r = 6-N^- 3AV - N v (3.2) 

for the model A and 

5 V = 6 - - 2Ntf - N v (3.3) 

for the Gribov model. 

In dynamical models, the 1-irreducible diagrams without the fields ijy vanish, and it 
is sufficient to consider the functions with ALt > 1- For the Gribov model, the functions 
built solely from the fields ip^ also vanish due to the symmetry (12.80 . (12.190 . and it is 
sufficient to take > 1 in (13.30 . For the model A, due to the reflection symmetry, 
it is sufficient to consider in (13.20 only the functions with even ALf + Afy,. With these 
restrictions, the analysis of the expressions (13.20 . (13.30 shows that in the both models, 
superficial UV divergences can be present in the following 1-irreducible functions: 

{ip ip) (5 = 2) with the counterterms ip'dtip, ^d 2 ^, ip 'if), 
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{ift'ipv) (5 = 1) with the counterterms ip^(vd)ip, ^(dv)ip. 

For the model A, superficial divergence can also be present in the function 

(ijji \jjipip) (5 = 0) with the counterterm ip^ip 3 , 
while for the Gribov model, divergences can also be present in: 

(ip'iftip) (5 = 0) with the counterterm ip ■, 

(ip^ip't/j) (5 = 0) with the counterterm (^) 2 ijj. 

The superficial divergence in the function (ip^ipvv) with 5 = and the counterterm 
i/j'i/jv 2 , allowed by the dimension for the both models, are in fact forbidden by the 
Galilean symmetry! All the remaining terms are present in the corresponding action 
functionals ( 12.131) and (12.141) . so that our models are multiplicatively renormalizable. 
It is important to stress that, since neither ip> nor if> are conserved in the full model, 
the both counterterms ip^(vd)i/j and i/j^(dv)i/j will be generated by the renormalization 
procedure. Thus to ensure the multiplicative renormalizability, one has to include the 
both such terms in the original action from the very beginning; hence the necessity to 
introduce the two terms with the velocity in ( I2.10p and ( 12. 17ft . 

The Galilean symmetry also requires that the counterterms ip^dtip and ip^(vd)ip 
enter the renormalized action only in the form of the Lagrangian derivative ip^ t ip, 
imposing no restriction on the Galilean invariant term ip^(dv)ijj. 

We thus conclude that the renormalized actions can be written in the forms 

for the model A and 

S«($) = S«($) + {Z^)\ - Z^i, 2 } (3.5) 

for the Gribov model, where the common part of the renormalized actions is 

S§{$) = ^ {-Z 1 V t + A (Z 2 d 2 - Z 3 r) - aZ 6 (d iVi )}^ + 5(v), (3.6) 

with 5(v) from (12TT6|) . 

Here A, r, g, u and a are renormalized analogs of the bare parameters (with the 
subscripts "o" ) and \x is the reference mass scale (additional arbitrary parameter of the 
renormalized theory). The renormalization constants Zi absorb the poles in e and £ and 
depend on the dimensionless parameters u, w, a and a. Expressions (13.4p ~ fl3.6p can be 

|| The arguments based on the Galilean symmetry are usually applicable to the velocity field governed 
by the Navier-Stokes equation, and generally become invalid for synthetic Gaussian velocity ensembles. 
It turns out, however, that for a Gaussian ensemble of the type (|2. 1 1|) with vanishing correlation time 
the Galilean symmetry of the counterterms indeed takes place. This issue, along with the consequences 
of the Galilean invariance for the renormalization, is discussed in |13j in detail. The proof given there 
is also applicable to the models (|2.13|) . (|2.14[) . 



(3.4) 
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reproduced by the multiplicative renormalization of the fields if) — > ipZ^, ^ — > ip^Z^ 
and the parameters: 

g = gn £/2 Z g , u = gfi e Z u , w = wfi c Z w , 

A = AZ A , t = tZ t , a = aZ T . (3.7) 

Since the last term S(v) given by ( 12 . 161) is not renormalized, the amplitude Do 
from ( I2.12p is expressed in renormalized parameters as Do = wqXq = wXfi^, while the 
parameters m and a are not renormalized: m = m, ao = a. Owing to the Galilean 
symmetry, the both terms in the covariant derivative Vt are renormalized with the same 
constant Z\, so that the velocity field is not renormalized, either. Hence the relations 

Z W Z\ = 1, Z m — Z a — Z v — 1. (3-8) 

Comparison of the expressions ( I2.13l) - (12.15l) and f]3.4j) ( f3T6|) gives the following relations 
between the renormalization constants Z\-Zq and (13. 7ft : 

Z\ = Z^Z^-t, Z 2 = ZiZ\, Z 3 = Z 2 Z T , Zq = ZiZ a (3-9) 
for the both models, 

Z^ = Z\Z 2 U Z 5 = Z\Z u Z^Z^t (3.10) 
for the model A and 

Z4 = ZgZxZ^Z^, Z5 = Z g Z\Z^Z^p (3-11) 

for the Gribov model. Resolving these relations with respect to the renormalization 
constants of the fields and parameters gives 

Z\ = Z 1 1 Z 2 , Z T = Z 2 1 Z 3 , Z a = Z x 1 Zq, Z u = Z 1 1 Z 2 2 Z^Z* > (3.12) 

for the both models (for the Gribov model, where we passed to the coupling constant 
u = g 2 with Z u = Z 2 g ), 

7 _ yl/2 yl/2 7- 1/2 rv\ _ yl/2 y—1/2 yl/2 /q 1 q\ 

^ip — A l A 2 ^4 i — A l A 2 ^4 ^O.lOj 

for the model A and 

Z„ = Z\' 2 Z^ 2 Zy\ 4 = Z^ZTZ- 112 (3.14) 

for the Gribov model. 

The renormalization constants can be found from the requirement that the Green 
functions of the renormalized models fl3.4p - fl3.6p . when expressed in renormalized 
variables, be UV finite (in our case, be finite at e — > 0, £ — > 0). The constants Z\- 
Zq are calculated directly from the diagrams, then the constants in (13. 7p are found from 
f l3.12p -f l3.14p . In order to find the full set of constants, it is sufficient to consider the 
1-irreducible Green functions which involve superficial divergences; all these functions 
are listed above. In the one- loop approximation, they are given in figured] for the model 
A and in figure [2] for the Gribov model. 

The solid lines with arrows denote the propagator (12. 7p . the arrow points to the 
field ip*. The solid lines without arrows correspond to the propagator ( 12.91) and the 
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Figure 1. One-loop approximation for the relevant 1-irreducible Green functions in 
the model (|2A3| . 



wavy lines denote the velocity propagator (vv)o specified in ( 12. lip . The external ends 
with incoming arrows correspond to the fields the ends without arrows correspond 
to ip. The triple vertex with one wavy line corresponds to the vertex factor ( 12.18)) . 

The minus signs in front of the bare terms for the functions ip) i-ir come from 
the Dyson equations. The first and the last diagrams for the functions (ijj'ijjijj)i-i r and 
(ip^ ip^ ip) i_ ir for the Gribov model have in fact two different forms, related by the mirror 
reflection, but they give equal contributions to the renormalization constants and are 
accounted by the factors of 2. What is more, the last diagrams for those functions 
involve closed circuits of retarded propagators (ipip^)o and therefore vanish identically. 
For the same reason, the last diagrams in the functions {ip^ipip)i-i r and for the model 
A and (il)'ipv)i-i r for the both models also vanish. 

All the diagrammatic elements should be expressed in renormalized variables using 
the relations f l3.7p -f pT9]) . In the one-loop approximation, the Z's in the bare terms 
should be taken in the first order in u = g 2 and w, while in the diagrams they should 
simply be replaced with unities, — > 1. Thus the passage to renormalized variables in 
the diagrams is achieved by the simple substitutions Ao — > A, r — > r, g — > gff^ 2 and 

Wo — > Wf/t. 

In practical calculations, we used the minimal subtraction (MS) scheme, in which 
the renormalization constants have the forms Zi = 1+ only singularities in e and £, with 



Effects of turbulent mixing on critical behaviour in the presence of compressibility 13 



l-ir 



{-iujZ 1 + Xp 2 Z 2 + XtZ 3 } + ^ 





^W)l-ir = ^ £/2 AZ 4 + 2 




-9V 



S/2 XZ, 




l-ir 



ipZi — iaqZ 6 





Figure 2. One-loop approximation for the relevant 1-irreducible Green functions in 
the model (f2A4)) . 



the coefficients depending on the completely dimensionless renormalized parameters u, 
w, a and a. The one-loop calculation is rather similar to that for the incompressible 
case (a = 0), discussed in [TJ] in detail (see also p2} [13]), and here we only give the 
results: 



w 



l--(3 + «) 



w 



1 a(a - I V 

, u (4a - 1) 
£ 4a 



Z 3 
3u 3w 



u 

1 + -, 

e 



1 H aa 



for the model ^4 and 
Zx = l 



u 



4e' 



1 



u w 



(3 + a) 



1 + 



1 + - - -a(a-l) 

^ s 

it (4a - 1) 

£ 8a 



1 + 



M W 



-aa 



(3.15) 



(3.16) 



for the Gribov model, with higher-order corrections in the couplings. To simplify the 
resulting expressions, we passed in (l3.15j) and ( I3.16P to the new parameters 

u — > m/167t 2 , w — > w/16ii 2 ; 
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here and below they are denoted by the same symbols u and w. 

For w — 0, the well-known one-loop results for the models (12.41) . (12.51) are 
reproduced (up to the notation). For the incompressible case (a = 0) expressions 
(I3.16P coincide with the results obtained earlier in |T4]. Furthermore, the first-order 
expressions (13.161) satisfy exact relations 

Zi(a) = Zi(l — a) for i — 1,2, 3; 

Z 4 (a) = Z 5 (l — a), Z 1 (a) — aZ 6 (a) = (1 — a)Z 6 (l — a), (3-17) 

which follow from the symmetries (12.81) . (12.191) of the Gribov model. These relations 
will play an important role in the analysis of the fixed points of the model (12.141) in the 
next sections. 

4. RG functions and RG equations 

Let us recall an elementary derivation of the RG equations; detailed exposition can 
be found in monographs [U E]. The RG equations are written for the renormalized 
Green functions G R = ($ • ■ ■ which differ from the original (unrenormalized) ones 
G = ($•••$) only by normalization (due to rescaling of the fields) and choice of 
parameters, and therefore can equally be used for analyzing the critical behaviour. 
The relation S R (&,e,fi) = iS($,e ) between the bare (I2.13p - ()2.15p and renormalized 
(I3.4p -( 13T61) action functionals results in the relations 

G(e ,...) = Z^z"?G R (e, /*,...)• (4.1) 

between the Green functions. Here, as usual, and ALt are the numbers of 
corresponding fields entering into G (we recall that in our models Z v = 1); e = 
{A , t , u , w , a , m , a } is the full set of bare parameters and e = {A, r, u, w, a, m, a} 
are their renormalized counterparts (we recall that «o = ct and mo = m); the dots stand 
for the other arguments (times/frequencies and coordinates/momenta). 

We use to denote the differential operation fid^ for fixed eo and operate on both 
sides of the equation (14.11) with it. This gives the basic RG differential equation: 

{V RG + + iY>7^t } G R (e, //,...) = 0, (4.2) 

where T>rg is the operation D M expressed in the renormalized variables: 

V RG = V^ + Mu + AA + Ma - lxV x - 1 T V T . (4.3) 

Here we have written T> x = xd x for any variable x, the anomalous dimensions 7 are 
defined as 

'jp = In Z F for any quantity F, (4.4) 

and the (3 functions for the dimensionless couplings u, w and a are 
U = V^u = u (-£ - 7„), 
(3 W = Vf.w = w (-£ - j w ), 

/3 a = V^a = -aj a , (4.5) 
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where the second equalities come from the definitions and the relations (13. 7p . The fourth 
/3 function 

(3 a = V^a = -aja (4.6) 

vanishes identically due to (14. 15j) and for this reason does not appear in the subsequent 
relations. 

The anomalous dimension corresponding to a given renormalization constant Z F is 
readily found from the relation 

If = (AA + (3 w d w + (3 a da) In Z F ~ - (eV u + £V W ) In Z F . (4.7) 

In the first relation, we used the definition ( 14.4=11 . expression ( 14. 3 p for the operation X> M 
in renormalized variables, and the fact that the Z's depend only on the completely 
dimensionless coupling constants u, w and a. In the second (approximate) relation, 
we only retained the leading-order terms in the /3 functions (14.5 j) . which is sufficient 
for the first-order approximation. The leading-order expressions ( 13.15ft . ( 13 . 1 6[) for the 
renormalization constants have the form 

%i W 

Z F = 1 + -A F (a,a) + — B F (a,a). (4.8) 
£ s 

Substituting (14. 8 j) into (14. 7p leads to the final UV finite expressions for the anomalous 
dimensions: 

jp = — uA F (a, a) — wB F (a, a) (4.9) 

for any constant Z F . This gives 

7i = 0, 72 = u>(3 + a)/4, 73 = -u, 74 = wa(a - l) 2 , 

75 = -3u + 3waa 2 , 7 6 = u{\ - 4a) /4a (4.10) 

for the model A and 

71 = -u/4, 72 = -u/8 + w(3 + a)/4, 73 = -u/2, 

74 = — u + wa(a — l) 2 , 75 = — u + waa 2 , 76 = u(l — 4a)/8a (4.11) 

for the Gribov model. 

The multiplicative relations ( 13. 12[) — ( 13. 14j) between the renormalization constants 
result in the linear relations between the corresponding anomalous dimensions: 

7a = 72 - 7i> 7r = 73 - 72, 

7a = 76-7i, lu = -7i - 2 72 + 74 + 75 (4.12) 
for the both models, 

27^ = 71 + 72 - 74, 27^ = 7! - 72 + 74 (4.13) 
for the model A and 

27v = 7i-74 + 75, 27J, = 7x + 74 - 75 (4. 14) 
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for the Gribov model. Along with (j4.10p . (14. lip , these relations give the final first-order 
explicit expressions for the anomalous dimensions of the fields and parameters. The 
exact relations (13.81) result in 



lw = -7A, lm = la = lv = 0, (4.15) 

while the relations (13.171) give 

7i(a) = 7,(1 - a) for i = 1,2,3; 

74(a) = 75(1 - a), 71(a) - aj 6 (a) = (1 - a)7 6 (l - a) (4.16) 
for the Gribov model. 



5. Attractors of the RG equations and scaling regimes 

It is well known that possible asymptotic regimes of a renormalizable field theoretic 
model is determined by the asymptotic behaviour of the system of ordinary differential 
equations for the so-called invariant (running) coupling constants 

V s gi{s,g) = Pi(g), gi(l,g) = gi, (5.1) 

where s = k//i, k is the momentum, g = {g{\ is the full set of coupling constants and 
g~i(s,g) are the corresponding invariant variables. As a rule, the IR (s — > 0) and UV 
(s — > 00) behaviour of such system is determined by fixed points g^. The coordinates 
of possible fixed points are found from the requirement that all the /3 functions vanish: 

PM = 0, (5.2) 
while the type of a given fixed point is determined by the matrix 

% = df3i/dgj\ g=g * : (5.3) 

for an IR attractive fixed points (which we are interested in here) the matrix Q is 
positive, that is, the real parts of all its eigenvalues are positive. In our models, the 
fixed points for the full set of couplings u, w, a, a should be determined by the equations 

fiu,w ,a,a(^*5 ^*) Cl*) 0, (5.4) 

with the (3 functions defined in the preceding section. However, in our models the 
attractors of the system (15. ip involve, in general, two-dimensional surfaces in the full 
four- dimensional space of couplings. First, the function (14.6p vanishes identically, so 
that the equation f3 a = gives no restriction on the parameter a. It is then convenient 
to consider the attractors of the system (15.11) in the three-dimensional space u, w, a; 
their coordinates, matrix (15. 3ft and the critical exponents will, in general, depend on the 
free parameter a. What is more, in this reduced space the attractors will be not only 
fixed points, but also lines of fixed points, which can be conveniently parametrized by 
the coupling a. Although the general pattern of the attractors appears rather similar 
for the both models, it is instructive to discuss them separately. 
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5.1. Scaling regimes for the Gribov model 

The one-loop expressions for the (3 functions in the model ( 13.5}) . (13 .61) are easily derived 
from the definitions (14. 5ft . relations (14. 12f) . (14. 141) and (14. 15f) . and explicit expressions 
iH]): 

f3 u = u [—6 + 3u/2 + w(3 + a)/2 — waf(a)] , 
(3 w = w [-£ + u/8 + w(3 + a) /A] , 

P a = u{2a - l)/8, (5.5) 

where the function f(a) = a 2 + (a — l) 2 achieves the minimum value /(1/2) = 1/2 at 
a = 1/2. For the set ( 15 .51) the equations (15 .41) have the following four solutions: 

(1) The line of Gaussian (free) fixed points: u* = = 0, a* arbitrary. 

(2) The point = 0, = 2e/3, a* = 1/2, corresponding to the pure Gribov 
model (turbulent advection is irrelevant). 

(3) The line of fixed points 

u* = 0, = 4£/(3 + a), a* arbitrary, (5.6) 

corresponding to the passively advected scalar without self-interaction. 

(4) The most nontrivial fixed point, corresponding to the new regime (universality 
class), both the advection and the self-interaction are relevant: 

4[ e (3 + n )-6f] 8[^4 + 3g] 

3(5 + 2a) ' 3(5 + 2a) ' 1 K ' 

We recall that a is treated as a free parameter, which the coordinates of the fixed points 
can depend on. For a = 1/2, the relations ( 14.121) and ( 14.161) give 7 a = 76 — 7i = 0, 
so that (3 a = —a^ a vanishes to all orders of the perturbation theory, irrespective of the 
values of the parameters u and w. It then follows that for the fixed points (2) and (4), 
the value a* = 1/2 is in fact exact and has no higher-order corrections in e and £. The 
expression for for the point (3), which corresponds to the exactly soluble Kraichnan 
model, is also exact, while the expressions for u* and for the fourth point are the 
first-order terms of infinite double expansions in e and £. 

Admissible fixed point must be IR attractive and satisfy the conditions > 0, 
W7* > 0, which follow from the physical meaning of these parameters. In the one-loop 
approximations, these two requirements in fact coincide for all the fixed points. In 
figure [3] we show the regions of IR stability of the fixed points (l)-(4) in the e-£ plane. 
For the points (l)-(3), the matrix (15.31) is either diagonal or triangular, and the stability 
is in fact determined by its diagonal elements fi, = dfii/dgi. 

For the point (1) one has 

so it is admissible for e < 0, £ < (region I in fig. [3]). Vanishing of the element Q a 
reflects the fact that the parameter a* for the point I is arbitrary, or, in other words, 
the point is degenerate. For the point (2) one finds 

Q u = e, Q w = -£ + m*/8 = -£ + e/12, Q a = u*/4 = e/6, 
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£=(3 + a)e/6 




£ 



Figure 3. Regions of admissibility of the fixed points in the model ()2.14|) . 



so it is admissible for e > 0, £ < e/12 (region II in fig. [3]). For the point (4) the 
elements dfi a /du, df3 a /dw vanish, so that the matrix Q appears block-triangular and 
the eigenvalue Q a = w*/4 decouples. The eigenvalues of the remaining 2x2 matrix for 
the couplings u, w look rather complicated, but the determinant is simple and equals 
w*u*(a + 5/2). This expression suggests that the matrix is positive for u* > 0, > 0, 
and careful analysis shows that this is indeed so. From the explicit expressions (15. 7\i 
we conclude that the point (4) is admissible for £ > e/12, £ < (3 + a)e/6 (region IV in 

fig. ED. 

The remaining sector III in fig. [3] is specified by the inequalities 

£>0, (3 + a)e-6£<0, (5.8) 

while the eigenvalues of the matrix (15. 3p at the point (3) are as follows 

4af 



Q u = -e + 2£ 



(3 + a 



f\Oj*)l £; ^0 







(5.9) 



with f(a) from (15. 5ft . The first inequality in (15. 8j) coincides with Q w > 0, while the 
second can be rewritten as 

4«e 



J) u (a* = 1/2) = -e + 2£ 



(3 + a 



/(1/2) > 0. 



(5.10) 



Comparison of equations (I5.9P and (15.10p shows that the condition Q u > can be 
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satisfied only in the region III, provided the coordinate a* satisfies the inequality 

(a* - 1/2) 2 < -L [6f - (3 + a)e] (5.11) 
oCtq 

(it is important here that /(1/2) = 1/2 is the minimum value of /). Thus we conclude 
that the third regime, for which the self-interaction is irrelevant, is admissible in the 
sector III, specified by the inequalities (15.81) . The corresponding attractor has the form 
of an interval on the line (15.61) . centered at a* = 1/2 and defined by the inequality 
(15. lip . This interval becomes infinite for a = and shrinks to the single point a* = 1/2 
for a — » oo. It should be noted, however, that for the purely transverse velocity field 
(a = 0) the term with the coefficient a® in (12.101) and all the subsequent expressions 
vanishes, and this parameter in fact disappears from the model. 

In the one-loop approximation (I5.5p . all the boundaries between the regions of 
admissibility are given by straight rays; there are neither gaps nor overlaps between the 
different regions. Due to higher-order corrections, the boundaries between the regions 
II and IV and between the regions III and IV can change and become curved. However, 
it can be argued that no gaps nor overlaps can appear between them to all orders. The 
similar situation was encountered earlier in the crossover between the long-range and 
short-range regimes of the Gribov model; see [21]. Here we only note that, in the present 
case, it is important that the models with u = or w = are "closed with respect to 
renormalization" so that the functions (3 U for w = and /3 W for u = coincide with 
the P functions of the Gribov and rapid-change models to all orders of the perturbation 
theory. 

It remains to note that the location of the boundary between the regions III and IV 
depends on the parameter a. For the a = 0, it is given by the ray £ = e/2, in agreement 
with the result, derived earlier in [H] for the incompressible case. When a increases, 
the ray rotates counter clockwise, and in the limit a — > oo it approaches the vertical 
ray e = 0, £ > 0. 

5.2. Scaling regimes for the model A 

The one-loop expressions for the (3 functions in the model (13.41) . (13.61) are easily derived 
from the definitions (I4.5p . relations (14.121) . (I4.13P and (I4.15p . and explicit expressions 
HO}: 

(3 U = u [—e + 3u + w(3 + a)/2 — waf(a)] , 
Pw = w K + w{3 + a)/ A] , 

P a = u(4a - l)/4, (5.12) 

where the function f(a) = 3a 2 + (a — l) 2 achieves the minimum value /(1/4) = 3/4 at 
a = 1/4. 

(1) The line of Gaussian (free) fixed points: = w* = 0, a* arbitrary. IR attractive 
for e < 0, f < 0. 

(2) The point u>* = 0, = e/3, a* = 1/4, corresponding to the pure model A 
(turbulent advection is irrelevant). This point is IR attractive for £ < 0, e < 0. 
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(3) The line of fixed points 

= 0, = 4£/(3 + a), a* arbitrary, (5.13) 

corresponding to the passively advected scalar without self-interaction. It involves the 
interval 

(a. - 1/4) 2 < -L [-e(3 + a) + £(6 - a)] , (5.14) 
loa£ 

which is IR attractive for £ > 0, e < £(6 — a)/(3 + a). This interval becomes infinite 
for a — > and tends to the finite value — (e + £)/16£ for a — > oo (note that the right 
hand side of the inequality f)5.14p is positive within the region of IR stability). 

(4) The fixed point, corresponding to the new universality class, where both the 
advection and the self-interaction are relevant: 

«,. = 4£/(3 + a), u. = 6{3+ $~J { *~ a \ a, = 1/4. (5.15) 

3(3 + a) 

This point is IR attractive for £ > 0, e > £(6 — oe)/(3 + a). 

When a increases, the boundary between the regions of stability of the regimes 
(3) and (4) rotates in the upper half-plane e-£ counter clockwise from the ray £ = e/2 
(a — > 0) to £ = — e (a — > oo). The regions of IR stability of the fixed points (l)-(4) in 
the e-£ plane are shown in figure |4]for some value of a < 6, when the boundary between 
the regions III and IV lies in the right upper quadrant; for a > 6 it moves to the left 
upper quadrant. 



6. Critical scaling and critical dimensions 

Existence of IR attractors of the RG equations implies existence of self-similar (scaling) 
behaviour of the Green functions in the IR range. In this critical scaling all the IR 
irrelevant parameters (A, \i and the coupling constants) are fixed and the IR relevant 
parameters (times/frequencies, coordinates/momenta, r and the fields) are dilated. In 
dynamical models, critical dimensions Ap of the IR relevant quantities F are given by 
the relations 

A F = d h F + KdF + lF, A UJ = 2- 1 *x, (6.1) 

with the normalization condition = 1; see e.g. [2] for more detail. Here 
are the canonical dimensions of F, given in table [TJ and 7^ is the value of the 
corresponding anomalous dimension f)4.4p at the fixed point: 7^ = 7jt(w*, to*, a#). This 
gives A^, = d/2 + 7^, A^t = d/2 + 7* t for the Gribov model, A^ = d/2 - 1 + 7^, 
A^t = d/2 + l + 7^,t f° r the model A and A w = 2 + 7* for the both models. Substituting 
the coordinates of the fixed points from sections 15.11 and 15.21 into the explicit one-loop 
expressions (14.101) . (14. lip for the anomalous dimensions and using the exact relations 
(I4.12p - fl4.14p . one obtains the leading-order expressions for the critical dimensions. They 
are summarized in tables [2] and [3] for the Gribov model and model A, respectively. 

The results for the Gaussian points (1) are exact; all the results for the fixed points 
(2) have corrections of order e 2 and higher. The result A w = 2 — £ for the fixed points 
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Figure 4. Regions of admissibility of the fixed points in the model (|2.13j) . The case 
a < 6 is shown; for a > 6 the boundary between III and IV lies in the left upper 
quadrant. 

Table 2. Critical dimensions of the fields and parameters in the model (|2.14l) . (|2.f 51) . 
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Table 3. Critical dimensions of the fields and parameters in the model (|2.13p . (|2.f 51) . 
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(3) and (4) is also exact, as follows from the general relations 7a = j w in (I4.15P and 
A u = 2 — 7^ in (16.11) . and the identity 7^ = £, which is a consequence of the fixed-point 
equation (3 W = with (3 W from (14.51) for any fixed point with 7^ 0. The other results 
for the fixed points (4) have higher-order corrections in e and £. 

The result A T = 2 — £ for the fixed points (3) is also exact: for u — 0, the 1- 
irreducible function (ip^ip) is given by the only one-loop diagram, which gives rise only 
to the counterterm if)^d 2 if}] cf. [17] for the pure Kraichnan's model. Hence the exact 
relations Z\ = Z% = 1 for u = 0. Then from the relations (13.91) and (14. 15ft it follows 
7 T = —7a = jw which gives 7* = — £ for the fixed points with w* 7^ 0. 

Fixed point (3) in table [2] illustrates the general fact that the dimensions A,/, 
and A^t for the Gribov model interchange under the transformation a — > 1 — a, as 
a consequence of the exact symmetry ( 12.81) . ( 12.191) . which results in the relations ( 13.171) 
for the renormalization constants. It remains to note that for any fixed point the critical 
dimension of the velocity field is given by the exact relation 2A,, = — £ + A w following 
from the form of its pair correlation function ( 12.111) . ( 12.121) . 

7. Discussion and conclusions 

We studied effects of turbulent mixing on the critical behaviour, with special attention 
paid to compressibility of the fluid. Two representative models of dynamical critical 
behaviour were considered: the model A, which describes relaxational dynamics of 
a non-conserved order parameter in an equilibrium critical system, and the strongly 
non-equilibrium Gribov model, which describes spreading processes in a reaction- 
diffusion system. The turbulent mixing was modelled by the Kazantsev-Kraichnan 
"rapid-change" ensemble: time-decorrelated Gaussian velocity field with the power-like 
spectrum oc k~ d ~^. The both stochastic problems can be reformulated as multiplicatively 
renormalizable field theoretic models, which allows one to apply the field theoretic 
renormalization group to the analysis of their IR behaviour. 

We showed that, depending on the relation between the spatial dimension d and 
the exponent £, the both models exhibit four different critical regimes, associated with 
four possible fixed points of the RG equations. Three fixed points correspond to known 
regimes: (1) Gaussian fixed point; (2) critical behaviour typical of the original model 
without mixing (that is, model A or Gribov model); (3) scalar field without self- 
interaction, passively advected by the flow (the nonlinearity in the order parameter 
in the original dynamical equations appears unimportant). The most interesting fourth 
point corresponds to a new type of critical behaviour (4), in which the nonlinearity and 
turbulent mixing are both relevant, and the critical exponents depend on d, £ and the 
compressibility parameter a. 

Practical calculations of the critical exponents and the regions of stability for all the 
regimes were performed in the one-loop approximation of the RG, which corresponds to 
the leading order of the double expansion in two parameters £ and e = A—d. It has shown 
that, for the both models, compressibility enhances the role of the nonlinear terms in the 
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dynamical equations. The region in the e— £ plane, where the new nontrivial regime (4) 
is stable (the corresponding fixed point is positive and IR attractive), is getting much 
wider as the degree of compressibility increases. In its turn, turbulent transfer becomes 
more efficient due to combined effects of the mixing and the nonlinear terms. 

Let us illustrate these general statements by the example of a cloud of particles, 
randomly walking in a nearly critical turbulent medium. The mean-square radius R(t) 
of a cloud of such particles, started from the origin at zero time, is related to the linear 
response function ( 12. 6p in the time-coordinate representation as follows: 

^)=/^ G ( t ,x), GM H«/(0,0)>, , = |4 (7.D 

For the response function, the scaling relations of the preceding section give the following 
IR asymptotic expression: 

G(t,x)=x- A -VF(^-,^-), (7.2) 

with some scaling function F. Substituting ( 17.21) into ( 17.1 ft gives the desired scaling 
expression for the radius: 

R\t) = ^-A,-A,t)/A. f (_r_j 

where the scaling function / is related to F in ( 17.2ft as 

f( z ) = J rfxi 2 - 2A * F(x,z). 

Directly at the critical point (assuming that the function / is finite at r = 0) one obtains 
from (17. 3ft the power law for the radius: 

R 2 (t) cxt n , Q = (d + 2 - - A^)/A„ = (2 - 7 ; - tW/A*; (7.4) 

the last equality following from the relations ( 16.11) . For the Gaussian fixed points the 
usual diffusion law R(t) oc t}l 2 is recovered. For the regimes (3) the exact result 
R(t) oc t 1 /^ 2 -?) is derived, which for the Kolmogorov value £ = 4/3 gives R(t) oc t 3//2 in 
agreement with Richardson's "4/3 law" olR 2 jolt oc i? 4//3 for a passively advected scalar 
impurity. For the other two fixed points the exponents in (17. 3p . (I7.4p are given by infinite 
series in e (points 2) or e and £ (points 4); the first-order approximations are readily 
obtained from the results given in tables [2] and [31 

For the case of incompressible fluid (a = 0), the most realistic values d = 2 or 3 
and £ = 4/3 lie in the region of stability of the passive scalar regime (3), so that 
the spreading of the cloud is determined completely by the turbulent transfer and is 
described by the power law ( 17.41) with the exact exponent = 2/(2 — £). As the 
degree of compressibility a increases, the boundary between the regions of stability of 
the regimes (3) and (4) in the e-£ plane moves such that the region of stability of the 
regime (4) is getting wider (see the discussion in section [5]). When a becomes large 
enough, these physical values of d and £ necessarily fall into the region of stability of 
the new regime (4), the advection and the nonlinearity become both important, and the 
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crossover in the critical behaviour occurs. The new exponent in ( I7.4p can be represented 
in the form 

Q^=Q^+5Q, 5Q = -( 1 ; + 1 ; t )/(2-0 (7.5) 

(we recall that A w = (2 — £) exactly for the both regimes (3) and (4)). For the Gribov 
model, substituting the one-loop expressions for the dimensions 7^ and 7*^ from the 
table [2] into (17.51) gives 

(3 + a)e - 6£ 

m = 3(5 + 2a)(2-0' (7 ' 6) 
Straightforward analysis of the expression (I7.6P shows that, within the region of stability 
of the regime (4), the quantity SQ is positive and grows monotonically with a. Thus the 
spreading of the cloud becomes faster in comparison with the pure turbulent transfer, 
due to combined effects of the mixing and the nonlinear terms, and accelerates as the 
degree of compressibility increases. For the model A, the exponents Q in ( 17. 4|) for the 
regimes (3) and (4) coincide in the one-loop approximation, but the main quantitative 
conclusion remains the same: compressibility enhances the role of the nonlinear terms 
and leads to the widening of the region of stability of the full-scale critical regime. 

Further investigation should take into account conservation of the order parameter, 
its feedback on the dynamics of the advecting velocity field (mode-mode coupling in the 
spirit of the model H of critical dynamics), non-Gaussian character and finite correlation 
time of the velocity statistics. This work remains for the future. 
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